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Transport properties in doped Manganite are studied by considering the simplified Double- 
exchange mechanism in the presence of diagonal disorder. It is modeled by a combination the 
Ising Double-exchange with Falicov-Kimball model. In this model, charge order ferromagnetic and 
transport quantities are presented by using Dynamical-mean field theory in the case of finite cou- 
pling Hund and disorder strenght. Those results suggest that, within this model the insulator phase 
can be stated in the range of below charge order temperature which is suitable with observation of 
experiment in manganites. 

PACS numbers: 71.27.+a 



I. INTRODUCTION 



An enormous interest in the family of the doped 
manganice oxides Ti_ I D I Mn03 ( in which T is trivalent 
rare earth and D is divalent alkaline) in perovskite 
structure have been recent renewed 4*2*&*As doping x 
and temperature T varied, the manganites show a very 
rich phase diagram involving phases with spin, charge 
and orbital orderiSiiiiln the cubic lattice, the five-fold 
degenerate 3d-orbitals of an isolated atom or ion are 
spilitted into a manifold of three lower energy (d xy , d yz , 
and d xz ), usually referred to as ti g one mixing with the 
surrounding oxygens is included, and two higher energy 
states (d 2 .2_ y 2 and d 3z 2_ r 2) called e g state. Electron 
t2 g are usually located but electron e g are itinerant so 
they create band energy. The correlation between ti g 



and e g electron is expressed by Double Exchange (DE) 
mechanism— inwhich two motions involving electrons 
moving from the oxygen atom to the Mn 4+ ions and 
other one from the Mn 3+ to the oxygen atom. The main 
feature of DE is cooperative effect where the motion 
of an itinerant electron favors the ferromagnetic (FM) 
ordering of local spin and, vice versa, the presence of 
the FM order facilitates the motions of the itinerant 
electron. The DE model has been studied in "one 
orbital" model because it is not only simplify but also 
contains enough degree of freedom to manifest the 
phenomena of DE mechanism which was envisioned in 
the early days of manganites studies. 
Recent experiment observed in the Manganites have 
shown that the DE model alone can not explain the 
quantitative features of the temperature dependence of 
the resistivity through the entire temperature range-. 
To improve the theoretical understanding in Manganites, 
there are two proposed resolutions of this problem. The 
first one is a large Jahn- Teller lattice distortion which 
causes a metal-insulator transition (MIT) via strong 
polaronic narrowing of the conduction electron band. 
However, the Jahn- Teller lattice distortion makes the 
charge-ordering phase occurs before the FM transition 
temperatureiwhich is not suitable with experimental 
observation at half-fillingASln the second proposed 



resolution, the insulating behavior is caused by a combi- 
nation of both magnetic disorder and non-magnetic ionic 
disorder. The magnetic disorder arises from the random 
DE factor cos(0/2) in the electronic hoping (6 is the 
angle between local spins) which call off-diagonal. The 
non-magnetic comes from the ionic doping of the D 2+ ion 
(i.e. from randomness at the chemical substituting T 3+ 
by D 2+ ) which leads to a random local potential for the 
charge carriers. This disorder always presents in doped 
Manganite which called diagonal-disorder that can lead 
to Anderson localization of the charge carriersiSTt is well 
known that, the diagonal-disorder can be modeled by 
the Falicov-Kimball (FK) modelifi Although the FK 
model is simple, it contains a rich variety of phases which 
are controlled by electron interaction^ii^Incorporating 
the diagonal-disorder of FK type into DE model, one 
may expect that a disorder-order phase transition could 
present. When phase transition occurs, a CO-FM 
phase may be stabilized at low temperature 1 ^ and 
the behaviors of transport coefficients can be more 
close with the experiment results. In this paper, the 
transport coefficients such as d.c. conductivity, thermal 
conductivity and thermal power of system are being 
studied by the Dynamical mean- field theory (DMFT) 14 . 
The DMFT have been extensively used for investigating 
strongly correlated electron system^It is based on the 
fact that, the self-energy depends only on frequency in 
the infinite dimension limit. Using DMFT we found 
that the system stabilizes inhomogeneous ferromagnetic 
state at low temperature and gives us a fulfill picture of 
behaviors of the transport coefficients. 
This paper is organized as follows: 



II. MODEL AND TRANSPORT COEFFICIENTS 

The system which we consider is described by the fol- 
lowing Hamiltonian 
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(1) correlation function at zero frequency, i.e., 



where c| (c i(T ) is the creation (annihilation) operator 
for an itinerant electron with spin a at lattice site i. 
The first term in the Hamiltonian (JIJ represents the 
hoping between the nearest neighbors sites, t is the 
hoping integral and is scaled with the spatial dimension 
d in the limit d — > oo as 14 , t — -^=. In the following 

we will take t* = 1 as the unit of energy^ 6*1 is the 
z component of the local spin at lattice site i, it takes 
values -§, -§ + 1, |, s? = | E CT <4cr z c iCT with cr z is 
the z- component of Pauli matrix is the spin of itinerant 
electron, ra/j is a classical variable that assumes the 
value 1(0) if site i is occupied by (is not occupied by) 
D 2+ ion. U is the disorder strength and is mapped 
onto the difference in the local potential which spilits 
energetical favor of Mn 3+ and Mn 4+ ion, the expectation 
value x = J2i( n fi) /N (N is the number od lattice sites) 
corresponds to the concentration of unfavorable Mn 4+ 
sites. The chemical potential fi controls the carrier 
dopings, while E nf controls the fraction of the sites 
having the additional local potential. The condition 
ri + x = 1 where n = ^2 ilT {nicr) /N is electron doping is 
used which determines Ef for each doping n. The first 
three terms in the Hamiltonian describe a simplified 
DE model which contains only the Ising-type interaction 
between the itinerant and local electrons. Within the 
DMFT the Ising type of the Hund interaction gives 
basically same results as the ones of the Hund interac- 
tion of clasical spins li^The simplication does not allow 
any spin-flip processes, which can be important at low 
temperature where spin-wave excitations may govern 
the thermodynamics of the system. However, within 
the DE model the DMFT calculations for quantum 
spins does not show significant differences from the 
one for clasical spinsii^iSiThe last two terms in the 
Hamiltonian Q describe a binary randomness of the 
T-site substitution. They together with the hoping term 
of itinerant electrons form the FK modeliZ^ 
Transport coefficients are calculated within a Kubo- 
Greenwood formalism- 9 , in which the dc-conductivity a, 
thermal power S and the thermal conductivity k can be 
expressed with the relevant correlation functions of the 
current operators. We have 19 
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where coefficients L n ,L 12 and L 22 are determined from 
the analytic continuation of the relevant current-current 
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Here L lJ (y) is the relevant current-current correlation 
functions 
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dre^ T (T r j(r)j(0)), (6) 



dre^(T T j(r)j Q (0)), (7) 



dre^ T (T T j Q (r)j Q (0)). (8) 



Where j and jq are the particle-current operator and the 
heat-current operator, respectively. The particle-current 
operator is defined by taking commutator of the Hamil- 
tonian with the polarization operator RifT-i^In our 
model, there are two inequivalent sublattice A and B so 
it is useful to divide Hamiltonian into four terms. 



H = H A + H l 



H 



AB 
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BA; 



(9) 



With H a = Uj2a n fa n aa ~ ^JY,a S a s a + E fJ2 n fa 

where a = A(B) for sublattice A(B) expresses the 
interactions of sublattice A(B) and H^b = ~iE,j4^! 
Hba — —t^2a ^CT a o- with aj., a a (b' a , b a ) are creation and 
annihilation operator for sublattice A(B) respond to the 
hopping between A and B site. After some algebraic 
calculations we obtain. 
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in momentum representation, where velocity 
u q = y q e(q) and e(q) is the dispersion of the 
non-interacting electrons. The heat-current jq can be 
determined by the energy current and the particle 
current j, jq = Je ~ /A]- The energy-current is defined by 
taking commutator of the Hamiltonian with the energy 
polarization operator ^\ Rj/ij (where H = J^i^i)- We 
obtain 



JQ = E v q[ e ( q ~^K <J VT + & qCT a qCT ) 

~\ E VF(q-q')(v q ,+v q )( aqCT 6 qV + 6 qCTaqV ),(ll) 
q.q'er 

where W(q - q') = ^ J2ii JcrS i ~ n f U)e^-^^. 
Using the equation of motion (EOM) technique^ we can 
write the heat-current operator in the form 
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Substubting j and jg from lfTU|) and (jTT|) into ©0 and 
! © , the transport coefficients are able to be calculated in 

ir> = lim - )v„ the canceling of the vertex correction at infinite dimen- 



q,(T 

x(aL(T)& qCT (r')+6+ (r)a qCT (r')). (12) 
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(13) 
(14) 
(15) 



where /(w) = l/(exp(w/T) + 1) is Fermi-Dirac distri- 
bution function and A acr (e,w) = — ilmG aCT (e, w) is the 
spectral function of the Green function for the intinerant 
electron for sublattice A or B. In that way, the transport 
coefficients are fully determined by the single particle 
spectral function. 



III. DMFT IN SDE WITHIN FK MODEL 

We solve model Q by the DMFT. The DMFT is based 
on the infinite dimension limit. In the infinite dimension 
limit the self energy is pure local and does not depend 
on momentum. By using the standard technique^ the 
matrix Green function of the itinerant electrons satisfies 
in the following form 



G" 1 (k,zc) = 



Wk + A»-££(w) -e(k) 

-e(k) w k + M-Sf(w) 



(16) 

where e(k) is the dispersion of free itinerant electrons 

on a hypercubic lattice, H^ B \u>) is the self energy for 
sublattice A(B). The self energy is determined by solving 
an effective single-site problem. The effective action for 
this problem is 

S a eft = - J ' dT J ' dT c t<r( T )Ga*( T ~ t ') c ^(t') 

-jdrj: [JS>* + M - Un fa] cUr)cUr) + E^f) 



where c aa — a ff (6 ff ) if a — A{B) and similar to c' aa , 
Quo (t — r ) is the Green function of a effective medium 
for each sublattice A or B. It plays as the base Green 
function of the effective problem. The local Green 
function also satisfies the Dyson equation. 



G-l{u) = g-l{u)-v aa {u). 



(18) 



The local Green function G aa (uj) of the effective 
single-site problem is solely determined by the partition 
function. It can be calculated by. 



G c 



SZ 



(lLU n 



aeff 



(19) 



where Z e g is the partition function and w„ = (2n + l)irT 
is the Matsubara frequency. Within the effective single- 
site problem, the partition function becomes 



*eff 



= Tr / Dc\, a Dc aa e 



(20) 



where the trace is taken over S z and n a f. This partition 
function can be calculated exactly. Indeed, the Dynamics 
of the local spin S z and impurity n a f involved in the ef- 
fective action i|17|) independently, we could take the trace 
over S z and n a f in partition function (|20|) independently. 
This is similar to DMFT solving the FK modelSi, we 
obtain 



exp I J2(HG**(^ n ) + Jnorn + fi — Un af ] - lniw n )e iw » e+ L e E f n «^, (21) 
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where m = — |, — § + 1, . . . | are projections of the local spin S z on z axit. 

Using Ea. l|19|) we obtain the local Green function on real axis after taking analytical continuation 

= £ -prr— — , (22) 

where 

Wma = ^ eff eX P j"^/ ^/M^^T^M^M + ^g^TO + /i - E/n. a /] 1 . (23) 



Equation (12211 is similar to the one with the classical local 
spin— It can be considered as a generalized formula of the 
local Green function for any value of the local spin. The 
self-consitent condition of the DMFT requires the local 
Green function of the original lattice. That mean, 

G aa {iS) = ^ > £ i G aa <&w). (24) 

k 

From equation we obtain 

G aa {uj) = (a a [ dep(e)- — -, (25) 

where a = A, B and a = B, A, ( aa = lu + fi — S"(oj). 
Eqs. (|18|H|22|I . and l|25|) form the complete set of equa- 
tions, which self-consistently determines the self-energy 
and the Green function for each sublattice A and B. 



IV. NUMERICAL RESULTS 

In this section we present the solution of the 
self-consistent DMFT Eqs. (0 , and fgft. The 
self-energy and the Green function are obtained by 
iterations. 21 With the obtained Green function, we 
calculate the d.c. resistivity, thermal conductivity, and 
thermal power by Eqs. l|T3 )l l(H )l . and 
In Fig. la we present the magnetization ^a(b) = 
YlieA(B)( n it ~ n ii)/N of sublattice A(B) as the func- 
tion of temperature and we also plot the dependence 
of charge-order parameter A = (jia — tlb) (where 
n A(B) = J2ieA(B)a n io'/^) 011 temperature in a case 
where the filling of the localized and intinerant electron 
is equal 0.5 with J — 8 and various of disorder strengths 
U. In this paper, we only concentrate in the case of 
normal manganites where the Hund coupling dominates 
with the disorder strenght that mean J# *C U. When 
U is large, it shows that, below critical temperature the 
magnetization of both-sublattices exists. They equal 
to each other until another critical temperature where 
the charge-order parameter exist. That mean, at low 
temperature, system is in the charge board CO-FM 
state. In this phase the charge-order coexists in the FM 
state which had been observedt^But with the small value 



of U, the CO-phase disappears, and only FM phase 
exists. With those coexistances of phases, the depen- 
dence of transport coefficients on temperature behaves 
as expected. In Fig. lb the behaviors of resistivity in 
temperature range is plotted. 

It shows that, in high temperature range, the depen- 
dence of resistivity on temperature shows the metallic 
behavior (dp/ dT > 0). Immediately below the transition 
temperature ferromagnetic Tc the resistivity declines 
with decreasing temperature due to the reduction of the 
magnetic contribution to electron scattering. Hence, 
at low temperature system is in FM phase. But when 
temperature decreases below CO transition temperature, 
the system changes the phase to insulator (dp/dT < 0) 
in the large U case. In the case of small U our system is 
in metallic through the entire temperature range which 
is the same results as in the SDE model—. 
The kinks at Tc and Tco temperature in the resistivity 
curve indicate the FM and CO transitions. The last 
two Fig.l(c and d) show the dependence of the thermal 
conductivity and thermal power vs temperature. They 
behaviors express as expected. At low temperature the 
thermal conductivity decreases and disappears when 
T = 0. When temperature increases, the thermal con- 
ductivity has a peak at the temperature which decreases 
with the decreasing of U, but with the large of T the 
thermal conductivity slowly increases and is reputed 
without depending on disorder strength U. The thermal 
power also behaves as expected. At low temperature, it 
is negative that because of the priority of hole-carrier 
particles and shows a sharp peak. When T is increased, 
the thermal power changes the sign and increases. Those 
results had been obtained in the FK model— 

Above results can be understood clearly if we plot 
the DOS of itinerant electron. In Fig. 2 the DOS for 
half filling nf = n e — 0.5 case with J = 8, U = 0.5 for 
various temperature is expressed. At low temperature 
(T < Tco) the DOS of each sublattice and each spin 
direction are distincted, that mean our system is in 
CO-FM phase. At temperature T co < T < T c the 
distinction between sublattice A and B disappears but 
ones for spin direction are still on distinction which 
expresses the FM phase only. When the temperature 
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FIG. 1: Temperature dependence of the magnetization mi(mj) for sublattice A(B) and charge-order parameter A(a), d.c. 
resistivity (b), thermal-conductivity (c), and thermal-power(d) in a case of half- filling (n e = n/ = 0.5) and J = 8.0 for different 
disorder strength U: U = 0.1 (solid line), U = 0.2 (dotted), U = 0.3 (dashed line), U = 0.4 (chain-dotted), and U = 0.5 (double 
chain-dotted). 




FIG. 2: Density of State of initerant electron in the case of of half-filling (n e = rif = 0.5), J = 8.0, and U = 0.5 for various 
temperatures 



is larger (T > Tc), all DOS of each sublattice and 
each spin direction are concided which dccribcd the 
Paramagnetic phase. In all temperatures above Too the 
Fermi-level states in band energy so our system is in 
metallic phase but with temperature is smaller than Tco 



the DOS are depleted at Fermi-level that mean there 
occurs the gap which expresses the insulator phase. 
The DOS for various U of spin up only, at low temper- 
ature are also expressed in Fig. 4. One can see that, for 
weak coupling the DOS is essentially unchanged from 
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FIG. 3: Temperature dependence of the magnetization m,i(mj) for sublattice A(B) and charge-order parameter A(a), d.c. 
resistivity (b), thermal-conductivity (c), and thermal-power(d) in a case of half-filling (n e +n/ = 1.0) U = 0.5, and J = 8.0 but 
for different concentrations: n = 0.5 (solid line), n — 0.6 (dotted), and n — 0.7 (dashed line). 
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FIG. 4: Density of State of itinerant spin-up electron only in 
the case of of half-filling (n e = rif = 0.5), J = 8.0, and some 
values of disorder strength: U — 0.1 (diamond), U = 0.4 
(square), and U = 0.5 (circle) (opened(filled) for sublattice 
A(B)) at low temperature T = 0.001. 
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FIG. 5: The critical temperature Tc (filled) and Tco 
(opened) as a function of disorder strength U for J = 8.0 
(square) and J = 6.0 (circle). 



the SDE model, but as U is increased further, the gap 
of DOS appears and its width increases. This result 
is precisely the kind of metal-insulator transition that 
Mott envisioned. 

The similar results for various itinerant concentration 
n e also are plotted in the Fig. 3. It shows that, when 
n e increases the Tco and Tc decrease which also is ob- 
tained in the infinite of Hund coupling Jii In Fig. 5 the 
dependence of the critical temperature of the FM and 
CO phase transition on disorder strength are presented. 
That means the disorder substantially decrease Tc of 
the FM transition. At the same time, with increasing 



U, critical temperature of the CO phase transition first 
increases, reaches its maximal value and decreases later. 
This behavior of critical temperature of the CO phase 
transition is similar to the one in FK modeli^. 
In the large strength of U the Tco and Tc are close in 
on the same value but Tco still is smaller than Tc ■ That 
means, the CO state is stabilized in the FM phase at low 
temperature. This results are observed in experiment 8 . 
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V. CONCLUSION 

In this paper, we have considered the transport 
properties in the SDE model with the combination 
of diagonal-disorder by employing the DMFT. It is 
found that the FM and CO states are stabilized at low 
temperature and expected behavior of the transport 
coefficients are expressed respectively. Those results 
improve the one for SDE model so they describe more 
suitable of the transport properties in the Manganites. 



Although, the manganites are too complicate a system 
to be completely described by this simple model, the 
obtained transport properties qualitatively reproduce 
the ones observed in the manganites. 
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